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ABSTRACT 

o 

We present angularly resolved spectra of T Tau North and South in the 3 fim 



(N 



water ice feature and K-band. Most of the water ice absorption lies along the line of 



sight toward T Tau South, confirming that it is viewed through stronger extinction. A 
decrease in ice-band absorption toward T Tau S between December 1998 and January 
2000, significant at the 2a level, was associated with an increase in its near infrared 
flux. Br7 emission is detected in T Tau North and South and H2 (2.12 fim) emission 
only toward T Tau South, consistent with previous studies of infrared companions to 
' T Tauri stars. Our results suggest that the near IR variability of T Tau S is probably 

caused by both variations in accretion rate and variable extinction along the line of 
sight. 
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S-j [ 1. Introduction 

T Tau is the eponymous pre-main sequence star (Joy 1945), but observational investigations 
over the last 20 years have revealed that it is an exceptional member of the class. Many of the 
young stars in the Taurus dark clouds are multiple systems, and in this respect T Tau is not atypical 
(Ghez et al. 1993; Leinert et al. 1993; Simon et al. 1995). It consists of T Tau North (T Tau 
N), the star seen in visible light at V ~10 (Herbig & Bell 1988), and T Tau South (T Tau S), 
discovered in the infrared by Dyck, Simon & Zuckerman (1982), at a projected separation ~100 
AU (0."7). T Tau S is itself a binary, with a projected separation of 7 AU (0."05) (Koresko 2000). 
Our observations do not resolve the 0."05 T Tau S binary, so for the purpose of this paper we will 
refer to it as T Tau S. 

The T Tau system is remarkable in part because the physical nature of T Tau S is unknown. 
It dominates the total system flux in the thermal IR, has a total luminosity two times that of T 
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Tau N (Ghez et al. 1991, hereafter G91), but has never been detected in visible light to limiting 
V magnitude of 19.6 (Stapelfeldt et al. 1998). T Tau S also shows pronounced variability in the 
infrared. G91 discovered that it flared by ~ 2 magnitudes in the near and thermal IR in 1990. 
They concluded that T Tau S experienced enhanced accretion, similar to an FU Orionis outburst. 
Kobayashi et al. (1994) and Simon et al. (1996) reported that the K band flux of T Tau S had 
decreased to the pre-flare value in the early 1990's. 

T Tau S is now regarded as a member of a small class of stars known as the infrared companions 
(IRCs)(Koresko et al. 1997, hereafter K97). IRCs are companions to visible T Tauri stars, are very 
luminous in the IR, and show pronounced variability. The near IR colors of IRCs imply a color 
temperature < 1500 K, lower than possible for a stellar photosphere. K97 suggest that the IRCs 
are young stars seen through strong extinction and are experiencing episodic accretion from their 
complex surrounding environment. 

Spectral diagnostics in the near infrared provide probes to study extinction and accretion in 
IRCs. Whittet et al. (1988) reported that the optical depth of the 3^m ice-band absorption feature, 
Ti ce , correlates well with extinction in visible light, A„. Sato et al. (1990) discovered the 3 fim ice 
feature in absorption toward the unresolved T Tau system. Leinert et al. (1996) monitored the 
ice-band optical depth toward the Haro 6-10 binary, an IRC system similar to T Tau, and found 
that both components varied in Tj ce , possibly attributable to variations in extinction along the line 
of sight. 

Herbst et al. (1995; 1996; 1997) detected H 2 emission in the v=l-0 S(l) line (2.12 fim) in the 
Haro 6-10 and UY Aur IRCs. The Herbst et al. (1996) observation of T Tau indicated peaks of 
unresolved emission at the positions of T Tau N and S. The absence of emission in the i>=2-l S(l) 
(2.25 /jm) transition in these sources suggests a line ratio consistent with a shock origin rather 
than fluorescence. The observed H 2 emission could be caused by accretion shocks from the infall 
of material onto the star. Prato & Simon (1997) used Br7 line emission as a surrogate for Ha to 
distinguish between the classical and weak-lined T Tauri stars and Muzzerole, Hartmann & Calvet 
(1998) showed that the Br7 line luminosity in classical T Tauris correlates well with accretion 
luminosity determined from blue continuum excess. 

To determine if the observed variability of T Tau S could be caused by variable extinction, we 
made angularly resolved observations of the near IR spectral energy distribution (SED) of T Tau 
N and S to test for variation in ice-band absorption on three occasions. Additionally, we obtained 
angularly resolved K-band spectra of T Tau N and S to study Br7 and H 2 accretion diagnostics, 
and if possible, to detect evidence of the photosphere of T Tau S. 

2. Observations 

Ice-band observations were taken on UT 1998 December 11, 1999 November 17, and 2000 Jan- 
uary 27 at NASA's Infrared Telescope Facility (IRTF) during scheduled observing time (December) 
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and as part of the IRTF service observing program (November & January). We used NSFCam, the 
facility near IR camera, which is equipped with a 256 x 256 InSb array. All of our observations 
were taken with the 0."056/pixel plate scale. To measure the spectral energy distribution (SED) in 
the ice-band, we used the K (2.2 fj,m) and V (3.8 jum) broadband filters and the circularly variable 
filter (CVF), which has a spectral band pass of 2%, set at 2.4, 2.9, 3.05 and 3.4 /j,m. In December 
1998 we also used a CVF setting at 3.6 /im and the H-band (1.65 jiva) filter. To freeze the seeing, we 
took 100 0.1 second exposures in the K and L' filters, 200 0.1 sec exposures at all of the CVF filters 
except 2.9 /im, The relatively weak flux of T Tau S in the H band and atmospheric absorption at 
2.9 /j,m required 400 0.1 sec exposures. 

Figure 1 is a contour plot of T Tau N and S at 3.05 /im, demonstrating the resolution of the 
system in the ice-band. 16 Tau and rj Tau were the photometric standards for the broadband and 
CVF observations, respectively. The zero-point fluxes for the CVF observations were determined 
using a polynomial fit to the broadband IR fluxes of Vega, as reported in Table 7.5 of Allen's 
Astrophysical Quantities (Tokunaga 1999). We determined the Point Spread Function (PSF) at each 
of the observed wavelengths by making shift-and-add (S/A) images from a series of short exposures 
of the photometric standards. The flux ratios of T Tau N and S at the observed wavelengths were 
derived using two different binary fitting procedures; the results were always consistent to within 
la. A discussion of the binary fitting routines will be included in a forthcoming paper (Beck et al., 
in preparation). 

Our K-band (2.05 to 2.4 (im) spectral observations of the T Tau system were made at the 
Univerisity of Hawaii 2.2 meter telescope on UT 1996 November 2 with KSPEC, the facility near 
IR spectrometer (Hodapp et al. 1994; 1996). This work closely followed S/A observations of T Tau 
in the K-band on UT 1996 October 26 at the IRTF which were subsequently used to flux-calibrate 
the KSPEC spectra. The KSPEC detector is a 1024 x 1024 pixel HgCdTe array with a plate 
scale of 0."167/pixel. We used the tip-tilt adaptive optics system with the telescope in its f/31 
configuration. With a slit width of 0."8, we obtained a resolution of ^760 at 2.2 fj,m. The slit was 
aligned at the position angle of the T Tau N-S system so that we acquired spectra of both stars 
simultaneously in a 35 minute exposure. 

We observed the spectrum of HR 1377, a B8 single star, immediately following T Tau to 
provide PSFs in the cross-dispersion direction as well as a reference spectrum to correct that of T 
Tau for telluric absorption and instrumental effects. The spectra of T Tau N and S were extracted 
by fitting, at each wavelength, model binaries constructed from the PSF of HR 1377, constrained by 
knowledge of the binary separation. The total amplitude and flux ratio of the binary model which 
yielded the best fit in the least squares sense was used to derive the individual spectral amplitudes 
of T Tau N and S at each wavelength. The spectra of T Tau N and S were divided by that of 
HR 1377 with the effects of its broad Br7 absorption removed following the procedure described 
by Greene & Meyer (1995). In order to restore the true form of the spectrum, we multiplied by a 
black-body curve of temperature 11,400 K, that of a B8 star (Drilling & Landolt 1999). 
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The flux ratio of T Tau S to N of these spectra at 2.2 fim, 0.18±0.01, is in excellent agreement 
with the value of 0.16±0.03 measured by the S/A observations in the K-band six days earlier at 
the IRTF. We applied the IRTF flux calibration to the spectra and show the final results in Figure 
2. Figure 3 is an expanded plot of the 2.10 to 2.20 /im region, which also shows a polynomial fit to 
the continuum. 



3. Results 

3.1. Absorption in the Ice-band 

Figure 4 presents the angularly resolved ice-band SED of T Tau N and S measured on UT 1998 
December 11. By fitting a smooth continuum distribution to the SED, we derive optical depths 
for the ice feature, Tj ce , of O.lzLO.l toward T Tau N and 0.5±0.1 toward T Tau S. If we consider 
the total unresolved system, we derive Tj ce = 0.3±0.1, consistent with Sato et al.'s (1990) value 
of 0.29±0.02. This agreement may, however, be fortuitous considering the photometric variability 
of the T Tau system and apparent variability of Tj ce (G91; Beck et al. 2000). Figure 5 shows the 
ice band SEDs of T Tau S for all three of our observations. The depth of the ice-band absorption 
decreased between December 1998 and January 2000 (Table 1). The observed change in n ce is 
significant at the 2a level and is associated with an increase in the near IR flux (Table 1). We do 
not detect a statistically significant change in ice-band optical depth toward T Tau N. 

It is useful to relate Tj ce and the more familiar extinction in the visible, A„. At this time 
relationships between Tj ce and A^ are available only for stars seen through molecular clouds (Whittet 
et al. 1988; Teixeira & Emerson 1999) and their applicability to the denser circumstellar regions 
considered here is uncertain (see also §4). Using Whittet et al.'s relation, Tj ce = [m^,, — Ay(0)] 
where m = 0.093±0.001 and A v (0) = 3.3±0.1, we determine A„ = 4±4 mag toward T Tau N 
and 9±2 toward T Tau S on 1998 December 11. On 1999 November 17 and 2000 January 27 the 
corresponding visual extinctions to T Tau S were 7±2 and 5±2 mag. We cannot use the relation 
derived by Teixeira & Emerson because it requires the full width at half maximum of the ice-band 
absorption and our CVF observations were too widely spaced to measure this. 

3.2. K-band Diagnostics 

The K-band spectrum of T Tau N (Figure 2) shows prominent Br7 emission and photospheric 
absorption lines of Na (2.208 /um), Ca (2.264 /xm), and the CO Av =2 bands. Visible light 
observations indicate that T Tau N has a spectral type of K0 and is obscured by A^ ~1 (Cohen & 
Kuhi 1979). J— H and H— K colors, derived from photometry of T Tau N appearing in Kenyon & 
Hartmann (1995), also indicate a value of A^ ~1. The photospheric absorption features detected 
in the K-band are consistent with a K0 spectral type, although the shallow depth of these features 
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in the spectrum of T Tau N requires a 2.2/zm continuum excess on the order of >70% of the 
stellar continuum (Figure 2). The slope of our T Tau N K-band spectrum implies an Ay of 10—12 
magnitudes for an early K to late G type stellar spectrum. Although complex nebulosity around T 
Tau is evident in several studies (e.g. Stapelfeldt et al. 1998), preferential scattering of the short 
wavelength light does not explain the discrepancy between the different estimates of A,,, since the 
near-IR colors also indicate low extinction. Rather, it is likely that spectral structure in the infrared 
excess emission across the K-band is resposible for the observed difference. These considerations 
will be discussed in more detail in Prato, Greene & Simon (in preparation). 

Figure 3 shows weak Br7 and H~2 v=l-0 S(l) emission lines toward T Tau S, at 3a and 2a 
significance, respectively, but no evidence for H2 emission in T Tau N. Table 2 gives the fluxes of 
these lines. Approximately three years after our spectra were obtained Davies et al. (2000) also 
observed the angularly resolved K-band spectra of T Tau N and S. Kasper(2000) and Kasper et al. 
(in preparation) provide a detailed discussion of these spectra. Their observation was made when 
T Tau S was in an elevated flux state (G91; Beck et al. 2000), with a 2.2 /im flux roughly triple 
the value we measured in October 1996. They detected Br7 emission toward both T Tau N and S. 
Their measured Br7 flux of T Tau S, 2.6±0.3 x 10~ 16 W m 2 s 1 , is ~ 2a higher than our value 
(Table 2) suggesting that its Br7 flux may vary. 

The spectrum of T Tau S presented in Davies et al. (2000) and Kasper (2000) shows no 
evidence of H2 line emission. This could indicate that its H2 line emission varied in the three years 
between our observations. The H2 emission could, however, arise from an extended circumstellar 
environment around T Tau S. In this case, the difference in line fluxes could be attributable to 
the different solid angles sampled, ~ (0.8") 2 for our observation, and (0.25") 2 for the Davies et al. 
spectrum. High resolution mapping observations are required to resolve this issue. 

4. Discussion 

The projected angular separation of TTau N and S is equivalent to only 100 AU at a distance of 
140 pc, but our observations show that there is significantly stronger ice-band absorption towards 
the IRC. G91 and Herbst et al. (1997) obtained the same result for absorption in the 10 /jm 
silicate feature and derived visual extinctions towards T Tau S of 4.6±2.5 and 5.4±0.5 mag for 
the two studies, respectively, consistent with our results in section 3.1. Our data from the three 
observations of the water-ice feature suggest an inverse correlation between the near IR flux and 
ice-band optical depth. While this has a natural explanation in terms of a decrease of obscuring 
material along the line of sight, there is a 1/6 probability that this ordering is the result of pairing 
of random data. Hence, further observations are required to determine if this correlation is real. 

If T Tau S is a young binary viewed through significant obscuration, the typical visual extinc- 
tions of 5 to 10 magnitudes derived from ice-band and silicate absorption features are not large 
enough to account for its observed near IR colors. K97 estimate that 35 magnitudes of extinction 
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are necessary to redden a 5200 K black-body to ~ 800 K, the characteristic near IR color temper- 
ature of T Tau S. However, 35 magnitudes of extinction should produce enormous ice-band and 
silicate absorption depths which are not observed. Two explanations seem possible. It is likely that 
the flux that is observed at the line center of the ice and silicate features does not originate from the 
stellar photosphere, hence the absorption does not trace all of the material that is present toward 
the system. Also, whether the size and composition of the dusty material near young stars, and 
the relations between and the ice and silicate optical depths, are the same as in the interstellar 
medium is unknown and is of considerable interest. 

Several possibilities for the geometry of the obscuring material in the T Tau system can account 
for all or a part of the extinction and its variability observed toward T Tau S. We may be observing T 
Tau S through a tenuous extension of the disk of T Tau N, which was detected in HCO + ( J = 3 — 2) 
by Hogerheijde et al. (1997), and in 3 mm continuum emission by Akeson et al. (1998). The 
extinction variability could be caused by inhomogeneous material in the disk of T Tau N moving 
through our line of sight, perhaps by Keplerian rotation. The T Tau N disk would have to lie within 
the 2:1 resonance to avoid truncation by tidal effects, corresponding to ~l/3 of the semi-major axis 
if T Tau N and S have equal masses (Artymowicz et al. 1991). This would imply a separation of 
T Tau N and S much greater than the projected separation of 100 AU. The obscuring material 
observed toward T Tau S may also be associated with a nearly edge-on disk around the 0."05 
binary, with rotation again responsible for the variability. Akeson et al. (1998) do not detect any 
mm continuum emission around T Tau S to a limiting flux of 9 mJy However, tenuous material 
that would escape detection at 3 mm could still be resposible for the extinction toward the T Tau 
S system. Another possible explanation for the observed extinction is the episodic embedding and 
accretion of T Tau S from a disk around the T Tau N - S system (K97). 

The observed extinction and variability toward T Tau S could arise from material which is not 
associated with a disk, but with the complex distribution of gas and dust surrounding the young 
triple system (Momose et al. 1996; Schuster et al. 1997; Weintraub et al. 1999). In this case, 
the absorbing material must have significant structure and large inhomogeneities in absorption on 
scales less than 100 AU to account for the present small extinction toward T Tau N. The irrregular 
dimming of T Tau N by 3-4 visual magnitudes in the early 1900's provides support for this scenario 
(Lozinskii 1949; Beck &: Simon, in preparation). 

There are indications that variations in extinction do not provide the entire explanation for the 
IR variability of T Tau S and the IRC phenomenon. The flare discovered by G91 was essentially 
gray, with a flux increase of ~ 2 magnitudes between 2 and 10 //m. Detection of Br7 emission 
associated with T Tau S in both our and the Kasper (2000) spectra indicates continuing accretion 
activity and evidence for variation in accretion rate. Additionally, variability of H2 line emission 
would indicate that the flows that shock molecular gas in the environment of T Tau S can vary on 
a timescale of a few years. 

The variability of ice-band optical depth and K-band emission features suggest that the near 
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IR flux variation observed toward T Tau S may be caused by changes in both extinction and 
accretion rate. Frequent, simultaneous and high precision measurements of the IR spectral energy 
distribution, including measurement of the ice and silicate features, and the Br7 and H2 emission 
lines, have the potential to clarify the nature of T Tau S and the IRC phenomenon. Improvements 
in image quality through the techniques of adaptive optics, and the powerful new IR cameras and 
spectrometers at the major observatories offer the means for rapid progress. 

We thank the telescope operators at the UH 88" and the IRTF, Chris Stewart, David Griep, 
Bill Golisch, Charlie Kaminski and Paul Fukumura-Sawada for their assistance with scheduled and 
service observations. We are grateful to T. Greene for providing the opportunity to acquire the 
T Tau K-band spectra. We thank R. White for helpful discussion, M. Kasper for the link to his 
Ph.D. thesis and information on his paper in preparation, and the referee for a helpful review. This 
research was supported by NSF Grants 94-17191 and 98-19694. 



-8- 



Table 1: Near IR Flux and Ice-Band Properties of T Tau S 



Date 


K 


L' 


T~ice 




(Jy) 


(Jy) 




1998 Dec 11 


1.07±0.08 


5.5±0.3 


0.5±0.1 


1999 Nov 17 


2.13±0.05 


6.0±0.1 


0.3±0.1 


2000 Jan 27 


2.1±0.2 


6.9±0.3 


0.2±0.1 



Table 2: Br7 and H2 Line Emission 



Species 


Wavelength 


T Tau N 


T Tau S 




(jim) 


(x 10~ 16 W m~ 2 s" 1 ) 


(x 10" 16 W m- 2 s" 1 ) 


H 2 


2.1218 


<1.5 (3cr) 


1.1±0.5 


H I 


2.1661 


7.5±0.5 


1.6±0.5 



-9- 



REFERENCES 

Akeson, R.L., Koerner, D.W. & Jensen, E.L.N. 1998, ApJ, 505, 358 

Artymowicz, P., Clarke, C. J., Lubow, S. H. & Pringle, J. E. 1991, ApJ, 370, L35 

Beck, T. L., Simon, M., Ghez, A. M., Prato, L. & Howell, R. R. 2000 The Near Infrared and 
Ice-band Variability of T Tau and Haro 6-10 in The Birth and Evolution of Binary Stars 
poster proceedings of the IAU Symposium No. 200, eds. Bo Reipurth and Hans Zinnecker, 
51 

Cohen, M. & Kuhi, L. V. 1979, ApJS, 41, 743 

Davies, R.I. et al. 2000, Adaptive Optical Systems Technology Conference Proceedings, SPIE 4007, 
ed. P. Wizinowich, 952 

Drilling, J. S. & Landolt, A. U. 1999, Allen's Astrophysical Quantities, A. N. Cox ed., (4 th edition: 
New York: Springer- Verlag) , 388 

Dyck, H. M., Simon, T. & Zuckerman, B. 1982, ApJ, 255, L103 

Ghez, A. M. et al. (G91) 1991, AJ, 102, 2066 

Ghez A. M., Neugebauer, G. k Matthews, K. 1993, AJ, 106, 2005 

Greene, T. P. & Meyer, M. R. 1995, ApJ, 395, 516 

Herbig, G. H. & Bell, K. R. 1988, Lick Obs. Bull., No. Ill (HBC) 

Herbst, T. M. Koresko, C. D. & Leinert, Ch. 1995, ApJ, 444, L93 

Herbst, T. M. et al. 1996, AJ, 111, 2403 

Herbst, T. M., Robberto, M. & Beckwith, S.V.W. 1997, AJ, 114, 744 
Hodapp, K.-W., Hora, J.L., Irwin, E. & Young, T. 1994, PASP, 106, 87 
Hodapp, K.-W., et al. 1996, New Astronomy, 1, 177 

Hogerheijde, M. R., Van Langevelde, H.J., Mundy, L.G., Blake, G.A. & Van Dishoeck, E.F. 1997, 
ApJ, 490, L99 

Joy, A.H. 1945, ApJ, 102, 168 

Kasper, M. E., 2000, Ph.D. Thesis, Max-Planck-Institut fur Astronomie 
Kenyon, S. J. & Hartmann, L. 1995, ApJS, 101,117 

Kobayashi, N., Nagata, T., Hodapp, K.-W. & Hora, J.L. 1994, PASJ, 46, L183 



-10- 



Koresko, C. D., Herbst, T. M. & Leinert, Ch. (K97) 1997, ApJ, 480, 741 
Koresko, C. D. 2000, ApJ, 531, L147 
Leinert, Ch. et al. 1993, A&A, 278, 129 

Leinert, Ch., Haas, M. & Weitzel, N. 1996, Disks and Outflows around Young Stars Conference 
Proceedings, CD-Rom Paper c411, ed. Beckwith, Staude, Quetz and Natta 

Lozinskii, A. M. 1949, Perem. Zvzd., 7, 76 

Momose, M., Ohashi, N., Kawabe, R., Hayashi, M. & Nakano, T. 1996, ApJ, 470, 1001 

Muzerolle, J., Hartmann, L. & Calvet, N. 1998, AJ, 116, 2965 

Prato, L. & Simon, M. 1997, ApJ, 474, 455 

Prato, L. 1998, Ph. D. Thesis, SUNY Stony Brook 

Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 618 

Sato, S., Nagata, T., Tanaka, M. & Yamamoto, T. 1990, ApJ, 359, 192 

Schuster, K. F., Harris, A. I. & Russell, A. P. G. 1997, A&A, 321, 568 

Simon, M. et al. 1995, ApJ, 443, 625 

Simon, M., Longmore, A. J., Shure, M. A. & Smillie, A. 1996, ApJ Letters, 456, L41 
Stapelfeldt, K. R. et al. 1998, ApJ, 508, 736 
Teixeira, T. C. & Emerson, J. P. 1999, A&A, 351, 292 

Tokunaga, A. T. 1999, Allen's Astrophysical Quantities, A. N. Cox ed., (4 th edition: New York: 
Springer- Verlag) , 150 

Weintraub, D.A. et al. 1999, ApJ 517, 819 

Whittet, D.C.B. et al. 1988, M.N.R.A.S., 233, 321 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



- 11 - 



Fig. 1. — A 3.05 /xm circular variable filter image of T Tau N and S from observations made on 
1998 December 10. The figure is oriented North up, and East to the left. The x-y coordinates are 
in units of 0."028, half pixels. The contours are 14, 28, 42, 56, 70 and 84% of the peak flux. The 
signal-to-noise is ~ 120 and the flux ratio (South/North) is 0.47±0.02. Overplotted in the upper 
right is a cross-cut of the image at an x position of 64. The apparent peaks at y=~25, 50, and 75 
are the first order Airy rings of the diffraction limited image. 

Fig. 2. — The angularly resolved K band spectra of T Tau N and S. The most prominent spectral 
features are identified. 

Fig. 3. — A magnified view of the 2.1 to 2.2 /j,m spectral region of T Tau N and S. H2 is detected in 
emission toward T Tau S, but not toward T Tau N. Fits to the spectral continua are overplotted. 

Fig. 4. — The total and angularly resolved near IR spectral energy distributions (SEDs) of T Tau 
N and S from data obtained on 1998 December 11. 



Fig. 5. — The observed change in the SED of T Tau S from 1998 December to 2000 January. The 
derived ice-band optical depths are listed in Table 1. 
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